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We consider a generalized Hubbard model with on-site 
and nearest-neighbour repulsions U and V respectively, and 
nearest-neighbour hopping for spin up (down) which depends 
on the total occupation rib of spin down (up) electrons on both 
sites involved. The hopping parameters are tAA, tAB and tsB 
for rib = 0,1,2 respectively. We briefly summarize results 
which support that the model exhibits s-wave superconduc- 
tivity for certain parameters and extend them by studying the 
Berry phases. Using a generalized Hartree-Fock(HF) BCS de- 
coupling of the two and three-body terms, we obtain that at 
half filling, for tAB < tAA =tsB and sufficiently small U and 
V the model leads to triplet p-wave superconductivity for a 
simple cubic lattice in any dimension. In one dimension, the 
resulting phase diagram is compared with that obtained nu- 
merically using two quantized Berry phases (topological num- 
bers) as order parameters. While this novel method supports 
the previous results, there are quantitative differences. 



I. INTRODUCTION 

Since the discovery of high T c much effort has been 
devoted to the study of the pairing mechanisms in highly 
correlated electronic systems. Two very well studied mi- 
croscopic models for such systems are the Hubbard and 
t — J models. While the search for superconductivity in 
the positive-fi/JIubbard model by numerical methods has 
failed so farocl, t — J-like models exhibit d-wave super- 
conductivity for certain parametersoE3. 

Although most of the high T c superconducting materi- 
als have a two-dimensional character, calculations in one 
dimension (ID) have been very useful. In one dimension 
the Hubbard and supersymmetric (J = 2t) t — -Ljaod- 
els can be exactly solved using the Bethe ansatzElrEj. It 
has been shown that for this particular case of the t — J 
model, the ground state consists of bound states with a 
gapless excitation spectrum, buJ-ibr J < 2t bound states 
exist for large enough densitiestj. 

Other strongly correlated models that introduce cor- 
related hopping interactions ace, subject of current re- 
search [14-51]. In 1989, Hirscho proposed a model for 
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the description of oxide superconductors by considering 
the holes in a nearly filled band as charge carriers. The 
Hamiltonian in standard notation reads: 

L 

f 1 - ^finj-o + (1) 

In contrast to the Hubbard model, the hopping am- 
plitude for single-particle hopping to a nearest-neighbour 
site depends on the occupancy of one of the sites involved 
in the process. 

A few years ago, a new integrable model of strongly 
correlated fermionic systems for hole superconductivity 
has been introducedli§E3. It is a modified version of 
Hirsch's model (Q) and its Hamiltonian reads: 

L 

j=l <J=±1 

(1 - Atn ]H1+a)/2 ,_ a ) (2) 

In ID, it has been solved using the Bethe .ansajtz in 
three parameter regimes: < At < 1, At < OliJiiJ and 
At > ilia. The two first regimes are related by a particle- 
hole transformation and present similar behaviour to the 
Hirsch model (|l|): there is a good quantitative agree- 
ment between the ground state energies, both models 
show gapless charge excitations and spin excitations with 
a finite gap for all densities and also dominating super- 
conducting correlations for densities larger than a critical 
valuaij. Instead, for the case with At > 1, both models 
(|l| and ||) behave quite differently: while the spin gap in 
the Hirsch model vanishes and there are no indications of 
superconductivity, the gap remains finite and supercon- 
ducting correlations are present for small doping in the 
Bariev model (Q). 

An important interaction that has not been consid- 
ered in the above models is the Coulomb repulsion. Wc 
consider a more complete Hamiltonian that includes on- 
site U and nearest-neighbour V interactions between the 
particles, and contains a more general kinetic term with 
different hopping parameters, depending on the total oc- 
cupation of the sites involved: 

H = - ®ij(4,- < 7 C 3-<r+ hc ) + U Yl ni 1 ni l + 

<ij>a i 

V ^2 n i n fi 
<ij> 
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<&ij = £aa(1 - n la )(l - Ujcr) + t BB n ia nj a + 

tAB(nia + n ]l7 - 2n icr n Jcr ) (3) 

The model (Q) corresponds to taking t AA + t BB — 
2t AB = 0, U = V = and At/2 = 1 - t AB jt AAl and 
the Hubbard model is obtained when V = 0, i^A = 
^bb = tAB- The case V = has been derived as 
the model that describes the— low energy excitations of 
intermediate-valence systemsEl and "hole" and cuprate 
superconductorsEirE3. The mode l with V was obtained 
for cuprate superconductocsEaEJ and proposed to de- 
scribe the benzene moleculecJ'El 

For t A s = V = CUand, \t AA \ = \t B B\ an exact solution 
of (||) exists in IDcZtEj, but this case is too peculiar, 
as discussed later. Hamiltonian (||) supplemented by a 
hopping of pairs of the type t p Y^Uj) c l ] c l 7 iQgJ, c j,i has 
also been solved in ID using Bethe ansatzEfEa for U = 
—2t p and t AA t BB = t AB - ^ s behaviour is similar to that 
of the Hubbard model. For negative U the model displays 
superconductivity. A more general model including this 
one and that of Bariev for particular parameters has also 
been solved using Bethe ansatzEJ. 

For t AB = 0,\t A A\ = \t BB \ = t,U,V ^ the 
phase diagram jaf the model (J3J) has been investigated 
previouslyE3'L3~E3 and some exact results have been found 
for the half- filled system (density of particles n = 1): 
a) in any arbitrary lattice in D dimensions with coor- 
dination number z = 2D the ground state is a Mott 
insulator (MI) with all sites singly occupied if U > z 
max(V, \tAA\ + \tBB\)', b) the ground state is a charge den- 
sity wave (CDW) ifV > max{U/z, U/2+\t AA \ + \t BB \)/2 
(simple cubic lattices) ; c) there is a region with mo- 
bile carriers which we call metallic (M) between these 
two phases, the M-MI boundary being at Um-mi = 
•z(|£a.a lilies |) (for D = 1 this state is a "non-conducting 
metal"E3M); d) for D = 1 the M-CDW boundary is at 
Vm-cdw = (U/2 + \t AA \ + \t BB \)/2. 

For < t AB < t AA = t BB = t, U, V ^ the phase 
diagram at n = 1 has been studied numerically in ID and 
2D and within mean-field approximationLlrE3. The high 
spin degeneracy of the MI phase for t AB = is lifted and 
gives place to a spin density wave (SDW) phase. There 
exists a metallic phase (the detailed nature of which will 
be discussed later) for small values of U and V, which 
shrinks as t AB — > t. When t AB = t (Hubbard limit), 
several results indicate that there is no M phase, except, 
eventually, on the second-order transition line between 
the CDW and SDW phases, that ends at a tricritical 
point. The position of this point is, to date, rWjUwell 
determined and should be around U ~ 2V ~ 4£E5~E3. 

The aim of the present work is to study, by numerical 
and analytical methods, possible existence of supercon- 
ductivity in the model (|3j) in regions of parameters for 
which exact results are not available. Special attention 
is paid to the occurrence of exotic p-wave triplet super- 
conductivity (TS). In particular, we show strong evidence 
that the extension to finite t AB of the above mentioned 



M phase displays TS. Some evidence of triplet supercon- 
ductivity exists in the extended Hubbard model (Eq. ^ 
with t AA =p-ip*iB — tBB) very near the line U = —2V for 
positive JTQET 

In Section 2 we briefly review the results which indicate 
that the model (Q) exhibits s-wave superconductivity (or 
dominant singlet superconducting correlations at large 
distances in ID) for certain parameters. At half filling 
we present evidence of a superconductor-insulator tran- 
sition. In Section 3 we explain the HF-BCS decoupling 
scheme and apply it to the electron-hole symmetric case 
t AB < t AA = t BB at half filling. We obtain a phase di- 
agram separating regions in which the stable phase is a 
CDW, SDW or triplet p-wave superconductor (TS). The 
latter is the stable one for sufficiently small U and V 
in a simple cubic lattice in any dimension. In Section 
4 we calculate the phase diagram in ID by. a numerical 
method recently introduced by two of ua!3, which uses 
topological quantum numbers as order parameters, and 
compare with the HF-BCS results. Section 5 contains 
the conclusions. 



II. S- WAVE SUPERCONDUCTIVITY 

The model (||) becomes simpler when the three-body 
part of it vanishes. The coefficient of terms of the form 

%tr C i,a n i-a n i-o IS H = ^AB - tAA - t BB . 

We concentrate first in the results for t% = 0. Us- 
ing a HF-BCS decoupling, it has been shown that the 
model leads to (extended) s-wave superconductivity for 
V = 0, small enough U and positive (negative)—^ = 
tAA — tAB for a more (less) than half-filled bandcil. In 
ID Japaridze and Miiller-Hartmann have calculated the 
correlation exponent K p using continuum limit theory 
and bosonizationH. They obtained that K p > 1 (su- 
perconducting correlations dominate at large distances) 
if 

U <U C = 8(t AB - t AA ) cos(?ro/2) - 6V (4) 

where n is the number of particles per site. 

The particular case ti = t AA corresponds to the rela- 
tion t AA =^rtfi B ,t AB — 0, and has been exactly solved 
for V — OEZrEI. Superconducting ^-paired states with 
off-diagonal long-range order (ODLRO) are part of the 
highly degenerate ground state. However, as a conse- 
quence of this degeneracy, the system does not display 
anomalous flux quantization (AFQ) for large ringsclJ. In 
other words, the Meissner effect is absent. 

Since, due to the rather pathological degeneracy of 
the ground state, the exact results do not allow to draw 
definite conclusions for t AB ^ while the bosonization 
results are expected to be valid only for weak interac- 
tions, we have carried out a detailed numerical study of 
the model for V = in rings of 10 and 12 sitesc3. For 
t AA = 1, t AB = 2, t BB = 3, 2/3 < n < 4/5 and small U, 
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there are clear indications of binding and AFQ (a ten- 
dency towards a periodicity of half a flux quantum in 
the energy as a function of flux). The calculation of the 
correlation exponent K p indicates that superconducting 
correlations dominate for U < U c , where U c ~ 9, 8 or 
6.5 for n — 0.4, 0.5 or 0.6 electrons per site respectively. 
These values are nearly 1.4 times larger than the cor- 
responding continuum limit results: 6.46, 5.66 and 4.70 
respectively, according to Eq. (|4|). For n = 0.8, tAA = 1, 
tAB = 1.5 and tss = 2 we obtain U c — 1.6, while Eq. (Q) 
gives U c = 1.24. 

The qualitative agreement between Eq. (Q) and our 
numerical results is lost at half filling. While the con- 
tinuum limit theory predicts that the system is an in- 
sulator for all positive values of [7e3, we find evidence 
of a superconductor-insulator transition as a function of 
U when t% — Hab — tAA — tsB = and any sign of 
ti = tAA — tAB (the sign can be-phanged using the sym- 
metry properties of the moddtlH). This agrees with, 
previous numerical work of A" p E3efl , 9-BCS calculation^ 
and recent results using slave bosonafj. For small posi- 
tive values of U we find in rings of 10 sites that the energy 
as a function of flux E{<f) has a form that suggests AFQ 
for tAA = 1, tAB = 0.5 and ^sp-i= 0. Furthermore, using 
topological quantum numbers!^ (as explained briefly in 
Section 4), we detect a transition from charge (p/ c ) and 
spin (7 S ) Berry phases (7 c ,7s) = (0,0) to (ir,n) (corre- 
sponding to the SDW phase) as U increases. We point 
out that the topological values (0, 0) are also obtained for 
the negative- U Hubbard model which displays singlet s- 
wave superconductivity. The value of U c separating both 
phases is small. For 10 sites we obtain U c — 0.075. From 
the size-dependence we estimate that this actually cor- 
responds to a lower bound. From the results of K p of 
Ref.E] and using symmetry arguments we estimate (re- 
member n = 1) U c = 0.5 for tAA — 1, tAB — 3/5 and 
tsB = 1/5 and U c = 0.3 for tAA — 1, tAB — 2/3 and 
tsB = 1/3, but careful finite-size scaling is necessary to 
give accurate values of U c . 

While the values of K p alone cannot distinguish be- 
tween singlet s-wave (even) and triplet p-wave (odd) 
superconducting states, the BCS results and the Berry 
phases (0,0) are indicative of the former. A demon- 
stration of the s-wave character was provided by the 
results of|-s.tochastic diagonalization by Michielsen and 
De RaedtcS which showed the presence of singlet-singlet 
quasi ODLRO in ID for tAA = 1, tAB = 1-4 and 
tBB = 1-8, n = 1.5 and U < 1, and also for tAA = tAB = 
tsB = 1, U = —4 and n = 1.5 (negative-?/ Hubbard 
model) . 

We discuss now the effect of the three-body term 
t 3 = 2t A B - tAA - t B B- When tAA = t B B the-piodel 
is electron-hole symmetric for bipartite latticesEj. For 
the electron-hole symmetric case with ta > 0, the mean- 
field HF-BCS decoupling of the two- and three-body 
terms leads to singlet s-wave (and also (i-wave in 2D) 
superconducting solutionscH. Howevepy,a positive t$ also 
favours CDW and SDW instabilities^. In 2D, the de- 



tailed (mean-field) calculation of the phase diagram for 
V = shows that the SDW is stable for positive U near 
half filling, while extended s-wave superconductivity is 
present in the ground state for small U or large doping. 
The rf-wave paramagnetic solution has always larger en- 
ergy than the other twoEJ. The electron-hole symmetric 
case with £3 < is analyzed in the next section. 

III. P-WAVE SUPERCONDUCTIVITY 

In this Section we discus the phase diagram of the 
model for the electron-hole symmetric case tAB < tAA — 
tBB and half filling, in simple cubic lattices in arbi- 
trary dimensions using HF-BCS. This decoupling leads 
naturally to p-wave superconductivity for small U and 
V. Although a real proof of the existence of this phase 
is lacking, further arguments given in the next Section 
support its existence in ID. For tAB = 0, exact re- 
sults have shown that for small U and V, the CDW 
and SDW phases become unstable because the system 
lowers its energy taking advantage of the kinetic energy 
terms tAA aad Jap, which are inactive in the CDW and 
SDW phasesEj'StJ. As mentioned in the Introduction, 
a phase diagram valid forjueveral lattices in arbitrary di- 
mensions was constructedEj. For tAB = V = , an exact 
solution exists in lDc2rE2l, but these exact results were 
unable to identify the nature of the third p hase for finite 
tAB ■ The numerical and mean-field results! 3 Ir 8 ! show the 
presence of mobile carriers and a non-vanishing Drude 
weight within the M phase in ID and 2D, suggesting 
that the system is a Luttinger liquid in ID and metallic 
in 2D for small U and V. However, a suggestion that this 
phase has dominant triplet superconductkig correlations 
at large distances was given only recentlyO. 

Note that if in the three-body part of the corre- 
lated hopping (see Eq. ||), the operator c\ a Cj^ a is re- 
placed by its expectation value r(assumed for this ar- 
gument independent of spin and nearest-neighbour pair, 
as in the non- interacting case), this term takes the form 
^ T Y,(%j) ni,- a nj,-a, with t 3 = 2t A B - t A A - t BB < 0. 
Thus, the three-body term reduces to an attraction 
of nearest-neighbour electrons with the same spin and 
triplet pairing is a natural consequence of it. 

The HF-BCS decoupling is the most convenient mean- 
field approximation to reduce the many-body terms of 
the Hamiltonian to one-body terms. It is a generaliza- 
tion of the procedure used by Foglio and Falicov for the 
normal caseEl A more direct way of obtaining the HF- 
BCS Hamiltonian (Hhfbcs) is to define the vacuum as 
the (unknown for the moment) Slater determinant that 
is the ground state of Hhfbcs- Then one should normal 
order the exact Hamiltopian with respect to this vacuum 
using Wick's theoremE3El Clearly the contractions that 
appear using this theorem are the HF-BCS expectation 
values. The exact Hamiltonian takes the form of the HF- 
BCS ground-state energy plus normal ordered one- and 
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many-body terms. Neglecting the latter, one obtains the 
HF-BCS Hamiltonian. 

Neglecting also for simplicity spin-flip expectation val- 
ues of the form (cj a Cj^ a ) (related to spiral spin struc- 
tures), we obtain Hhfbcs using the following approxi- 
mations for the two- and three-body terms of H: 



+(cI T cJi( c ijCiT> + h.c.) - |(c a c lT )| 2 , 



(•5) 



symmetry breaking of the expectation values entering 
Eqs. (5-8) with respect to the unperturbed system. Since 
singlet s- and d-wave solutions do not exist for tAB < 0c3, 
we have looked for triplet p-wave superconductivity (TS). 
Based on symmetry properties expected for the ground 
state in ID (explained at the end of the next section), 
we assumed that only the S x = ±1 components of 
the triplet order parameter do not vanish like in the 
Anderson-Brinkman-Model phase of superfluid 3 _ffeEl 
For this case we have (c||C^) for all i, j in one lattice 
direction and a vector S connecting nearest-neighbours 
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) = 0. As a consequence also 
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-(4t c jt)(( n u) + ( n Ji)) 

-(4Ai)( c n c 3i) ~ < c I T c a)< c a c jT>' 

and the same interchanging spin up and down. Choosing 
one spin orientation for the sake of clarity the three-body 
terms are replaced as: 

* cWWXr^) + 1(4, c^)! 2 - |(4 iSi )| 2 ) 

t J 



(nil))) 



- c Ii c il ( < c i T c iT > < c ii c iX > + ( c h c jl)( c il c tf)) 

~ c ]i c n(( c h c 3i){4i c ji) + ( c h4i){c n c n )) 
+mi((4i c jT)( n ii) + (4i4i)( c n c 3t)) 

+njl((4i c n)( n *l) + (4Ai)^i c J^ 

+4Ai^ C il C 3l)( n 3l) ~ ( C jl C 3l)(4jl C il)) 

+4t c ji(( n il)( c il c jl) - < c U c iT>< c IiSi» 
+c ll c n ((c^cl l ){n n ) - {c\A n ){4i C 3i)) 

+ c Ji c jT« c iT c 5i>< n ii> - (4Ai)(4i c n)} 
+(4i4i( c 3i c ii) + c 3i c iMi4i>)^% c 3i) 

+2[(4^3i)(\(4i c 3i)\ 2 - (miHnn)) 
+(4r4i)(( c ji c 3f)(4i c ii) - ( c u c jr)( n 3i)) 
+(4Ai)(( c n c 3i)(4i c ji) - ( c 3i c rt)( n n)) 

-(4r^)\(4A)\ 2 i 



Each thermodynamic phase of the model in the 
HF-BCS approximation is characterized by a different 



depends on direction. We also considered the usual 
CDW and SDW phases for which the cubic lattice in 
D dimensions is divided into two equal interpenetrating 
|(c io .c- ,}| 2 ]}, (6) sublattices A and B, in such a way that the nearest- 
neighbours of any site of A lie in B. In this case, for 
the SDW (rii, ff ) = (1 + mcre 4Q R 0/2 while for the CDW 
{n i: a) = (l + Ae iCJ Ri )/2 with Q = (tt,w, ...,%) and m and 
A order parameters. The resulting one-particle Hamilto- 
nian has the form of the non- interacting one with a renor- 
malized effective hopping t e ff, plus a symmetry break- 
ing perturbation which, at half filling, depends on two 
or three parameters to be determined selfconsistently. 
These parameters are the corresponding order parameter 

^\ip, m or A) and the different values of (4+5 a Ci ^) : one 
for CDW, SDW, or TS in ID, two for TS in' more than 
ID (for TS in more than ID t e ff becomes anisotropic). 

The phase diagram has a natural energy scale which 
we call E: the absolute value of the energy of the half- 
filled non- interacting case for hopping \t$\ = tAA + tBB — 
2tAB, counting both spins. The ratio E/\ta\ depends 
only on the dimension of the simple cubic lattice and is. 
4/tt = 1.273 in ID, 16/tt 2 = 1.621 in 2D, 2.005 in 3DEll 
and 0.798 in ooDEil. The only point in the U — V phase 
diagram for which the paramagnetic solution is stable is 
the triple point (Ut,V t ) = (E,E/z) where z = 2D is the 
coordination number. If and only if V < Vt, for any 
U, the paramagnetic solution is unstable against the TS 
phase; if and only if U > Ut = E, it is unstable against 
SDW; and if and only if V > {E + U) /2z it is unstable 
against the CDW. These boundaries allow us to establish 
that, in any dimension, inside the region bounded by 
the dashed line in Fig. 1, the TS is the stable phase. 
Also the CDW-TS boundary lies between the dashed and 
dot-dashed lines of Fig. 1. With an adequate change 
of variables, it can be shown that the self-consistency 
problems for the CDW and SDW take the same form on 
the line U = zV . Thus, in any dimension, for U > E, 
U = zV is the boundary of the CDW-SDW first-order 
transition, independently of the other parameters of the 
model. 

The above results are quite general and valid in any 
dimension. The region of stability of the TS extends al- 
ways beyond the dashed lines of Fig. 1 and the precise 
location of the boundaries of the TS depends on tAB and 
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the dimension. In Fig. 1, the specific case of ID and 
tAA — tBB — 1, tAB — 0.2 is shown. The change in slope 
of the SDW-TS boundary is due to a metamagnetic first- 
order transition inside the SDW from large to small order 
parameter as V increases. A similar situation occurs in 
the CDW-TS boundary. Near the triple point, and par- 
ticularly for tAB ~ 1) all order parameters become very 
small, and the numerical method used to solve the self- 
consistency equations breaks down. Similarly, in ID, for 
V = and tAB — 0.097, the effective hopping vanishes in 
the paramagnetic phase and the HF-BCS approximation 
becomes invalid for tAB ~ 0.1 or smaller. 

One expects that the HF-BCS results are reliable for 
small values of the interactions (small U, V and 1 — tAB, 
with tAA = tBB = 1) and D> 1, supporting the existence 
of the p-wave TS for these parameters. In ID there is no 
true LRO and the HF-BCS results only have qualitative 
validity. A comparison with more accurate methods in 
ID is made in the next section. 



IV. PHASE DIAGRAM OBTAINED FROM 
TOPOLOGICAL TRANSITIONS 

In recent years, the concept of Berry phase was a sub- 
ject of great interest in a variety of fields in physics. Zak 
has shown that it can be used for labeling energy band in 
solidsEa, and subsequent work showed that changes of po- 
larization are pranortional to the corresponding change 
in a Berry phaseEa. Ortiz and Martin have-generalized 
these concepts to a many-body ground stated, and Resta 
and SoreU* used this concept to identify a ferroelectric 
transitions. This many-body Berry phase is simply the 
phase captured by the ground state in a ring of L sites 
as the boundary conditions c\ +L a — e^c^^, complete a 
cycle from a flux <j) = to <j) = 27i@fiS We call this 
phase the charge Berry phase 7 C . 

While previous work assumed always that the ground 
state is non-degenerate for all (j> (except at isolated points 
of the parameter space at which j c is indefinite) , two of 
us have recently generalized this concept for the case in 
which there is a-crossing of levels in the ground state as 
a function of This is the case of the superconduct- 
ing phases, for which there is AFQ as a consequence of 
the crossing energy levels. Furthermore, a "spin" Berry 
phase 7 S was defined as that captured by the ground state 
in the cycle < <j) < 2tt varying the boundary conditions 



as C i+L a = elCT ^ c ] o- with cr=l (-1) for spin up (down)Ej. 
Due to the inversion symmetry of the model Eq. (|^) . the 
Berry phases are quantized and in ID can take only the 
values and ir (modulo 2ir). Thus, it turns out that 
7 c /7r and Js/tt are topological quantum numbers which 
are related to changes in the total polarization and the 
difference between polarizations for spin up and down 
respectively^. ._. 

Defining a vector 7 = (7 c ,7s) 5 it is easy to seec3 that 
for a CDW (SDW) with maximum order parameter one 



has 7 = (0,0) (7 = (7r,7r)), and these values are con- 
sistent with the change in up and down polarizations 
when the electrons of spin down of one phase are moved 
to the other sublattice to form the other phase. Since 
7 jumps in ir only at the phase transitions, the values 
7(CDW)=(0,0) and 7(SDW)=(tt, tt) are valid for any 
non-zero magnitude of the corresponding order parame- 
ters. These values are also consistent with the canonical 
transformation. 

4T = c i,T,4; = (-l) i 4 (9) 

under which the CDW and SDW are interchanged and 
the Berry phases modulo tt transform as: 



7 = 7 + 7T 



(10) 



Since the ground state for U = V = is invariant 
under Eq. (Q), Eq. ( [To| ) implies that for the third phase 
present in the diagram 7 C (T5) = j s (TS) + n (numeri- 
cally it turns out that j(TS) — (0, n j), and thus, at least 
one of the topological numbers 7 c /7r, 7 s /7r jumps at each 
boundary. Thus, this method combined with finite-size 
scaling is able to determine accurately the phase diagram 
(see Fig. 2). 

The best previous numerical methods to determine 
phase diagrams of this type were based on the size de- 
pendence of SDW and CDW order parameters, ant},!-! 
probability densities in a Monte Carlo sampling 
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However, these quantities as well as different correlation 
functions vary smoothly at the transijtipn-and it is very 
difficult to obtain accurate boundaries£j~L3. Instead, the 
use of topological numbers as order parameters necessar- 
ily leads to sharp transitions. 

In Fig. 2 we show the phase diagram of the model 
Eq. (|3|) in ID obtained with the above mentioned method, 
for two values of tAB, and compare it with the corre- 
sponding— HF BCS results and with the exact result for 
t ^B —> 0E3. Although, as mentioned in the previous sec- 
tion, the HF BCS results are not expected to have quan- 
titative validity in ID for large values of the interaction, 
for small tAB the resulting phase diagram is in reasonable 
agreement with that obtained with the above explained 
more reliable numerical method. The results of the lat- 
ter tend to the exact phase diagram as tAB ~ * . The 
comparison also shows that the TS phase in ID extends 
beyond the expectations of the HF BCS results. This is 
particularly clear for tAB ~ tAA — tBB and U ~ 2V, 
and is probably related to the fact that, in the contin- 
uum limit theory, the backscattering and Umkiapp terms 
coming from the U and V terms of the Hamiltonian (ul- 
timately responsible of the insulating-behaviour), nearly 
cancel each other on the line U = 2VoE3. 

We should note that the method of the topological 
transitions alone, is not able to identify the nature of the 
TS phase, since different phases can have the same topo- 
logical numbers. For example, we find 7 = (0, 0) not only 
for the CDW, but also in the negative-?/ Hubbard model 
with small negative V, for which s-wave singlet pairing 
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occursEHlEZI. We also find it for our model (Eq. ([§) with 
ta = V = and small U, for which also singlet super- 
conducting correlations dominate at large distances, as 
explained in Section 2. 

In the following we summarize the evidence in favour 
of the triplet superconducting phase in ID. 

a) Superconductivity: 

• Since each time the Berry phase is evaluated, the 
ground state energy as a function of the flux E(<p) 
is computed at the same time, we have checked that 
there is a tendency to AFQ (i.e. E{4>) ~ E + 
inside almost all the TS phase. 

• Our previous calculations of the correlation expo- 
nent K p for V — show the dominance of super- 
conducting correlations at large distances (K p > 1) 
for U < U„, and an opening of a charge gap for 
U > C/ C t3u3. The value of U c has been estimated 
as U c ~ 3.5 for t AB = 0.2 and t AA = t B B = 1^3 in 
good agreement with the result shown in Fig. 2. Ear 
t AB = 0.6, from K p one obtains U c = 2.05±0.05E3, 
while from, the topological transition we obtain 

u c = 2.\m. 

• For a "metallic" gapless phase with K p < 1 (as 
in the non-interacting -ease), one expects unde- 
termined Berry phaseal3'E3 instead of the result 
7 = (0, tt) for the TS phase. 

b) Triplet character: 

• The HF BCS decoupling in any dimension leads to 
unstable singlet and stable triplet superconductiv- 
ity. 

• For known cases of singlet superconductivity (men- 
tioned above) we obtain 7 = (0, 0) in contrast to 
the result j(TS) = (0,tt). 

• We find that in rings of 10 sites the nearest- 
neighbour triplet-triplet correlations are larger 
than the singlet-singlet ones at the largest distance 
in the ring (5 sites). 

• Numerically, we obtain a non-degenerate ground 
state inside the TS phase which should transform 
into itself under the symmetry transformation (|^) 
for U = V = 0, since the Hamiltonian is invariant 
at that point. An ordinary BCS singlet solution 
transforms into a SDW under this transformation, 
leading necessarily to a degenerate ground state for 
U = V = 0. Instead, using c'\ , = C-k+Q,i, |0) = 

life C 'I 1 an d cos(fc Q + Q a ) = — cos(fc Q ), it can be 
easily checked that our p-wave triplet BCS solution 
in the ^-direction (see Section 3 after Eq. (8)): 

\ Px >= n K+^ c tAv,)i°> (ii) 

fcfT,fca; >0 

is invariant under Eq. (^). 



• Finite-size scaling is consistent with the absence of 
a spin gap inside the TS phase. Continuum-limit 
theory predicts that triplet superconducting corre- 
lations dominate at large distances (due to loga- 
rithmic corrections) when K p > 1 and a spin gap 
is absentE3. 



V. SUMMARY AND DISCUSSION 

We have studied the occurrence of superconductiv- 
ity in the generalized Hubbard model Eq. (||). This 
model contains the most important physical ingredients 
expected to describe transition metals in generalc3 or 
cuprate superconductorsEiTcJ. In ID, most of the pa- 
rameters of the model appropriate for .systems like trans- 
polyacetylene have been estimatedEjO and are consis- 
tent with several regions of parameters for which we find 
superconductivity. 

When the three-body term £3 = 2t A B — t AA — tsB 
vanishes, singlet s-wave superconductivity is expected 
for electron densities per site n ^ 1 or very small val- 
ues of U. For n not too near one and V = 0, Eq. (^) 
describes qualitatively the values of U below which su- 
perconducting correlations dominate at large distances in 
ID. However, numerical results suggest that Eq. (Q) un- 
derestimates these values by a factor ~ 1.4. For n = 1, 
the region of superconducting behaviour is very small. 

Instead, for n = 1 and t A B < t AA = tsB = 1> we 
find evidence of a p-wave triplet superconducting phase 
(TS) for small values of U and V. The Hartree-Fock BCS 
approximation gives a stable TS for sufficiently small U 
and V in any dimension. In ID, the numerical method 
of the Berry phases, described in the previous section, 
predicts a stable TS even for small 1 — t A B an d large 
values of U ~ 3 near the line U — 2V. 

While a definite proof of the triplet character of this 
phase does not exist so far, there are several arguments 
in favour if it enumerated at the end of the previous sec- 
tion. In contrast to the previous case of s-wave super- 
conductivity (2t A B = t AA + tsB) mentioned above, for 
tAB < t AA = tsB, the effect of doping-aeems to weaken 
the TS in favour of the SDW at V = OB 

Experimentally, evidence of a triplet odd-parity su- 
perconducting order parameter exists in some quasi- 
1D organic conductors^ In ID, our solution is 
compatible with a nodeless superconducting gap, as 
evidenced by thermal conductivity measurements in 

(TMTSF) 2 CIC>M 

The effect of phonons, particularly in ID, can stabilize 
insulating CDW states, or favour s-wave superconducting 
states. However, based o&,the general study performed 
in ID using bosonizationcH, we believe that our results 
remain qualitatively valid in presence of a small or mod- 
erate electron-phonon coupling. In addition, if phonons 
are treated in the antiadiabatic approximation, the model 
Eq. (3) retains its form [21b]. 
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Figure 1: Mean- field (HF BCS) phase diagram of the 
model (||) for one particle per site and tAA = tBB> tAB- 
The full line is the CDW-SDW boundary in any dimen- 
sion. Inside the region bounded by the dashed lines, 
the triplet superconducting (TS) phase is the stable one. 
The CDW-TS boundary lies between the dashed and 
dot-dashed lines. E/2(tAA ~ tAB) is a number of or- 
der one which depends on dimension (see text). Solid 
circles (squares) are points on the CDW-TS (SDW-TS) 
boundary in ID for tAA = ^bb = 1 and tAB = 0.2. 

Figure 2: Phase diagram of the model (g) in ID for 
one particle per site, tAA = tsB = 1 and two different 
values of tAB- Solid (open) circles are determined from 
the jump of 7 C (7 S ). Solid squares, dotted and dot-dashed 
lines are HF BCS results. The dashed lines are exact 
results in the limit tAB — > OEj. Full lines are guides to 
the eye. The small difference between the results from 
both Berry phases at the CDW-SDW boundary is due to 
finite-size effects. 
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